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ABSTRACT 

We present empirical relations for determining the amount by which the effective tempera- 
tures and radii — and therefore the estimated masses — of low-mass stars and brown dwarfs are 
altered due to chromospheric activity. We base our relations on a large set of low-mass stars in 
the field with Ha activity measurements, and on a set of low-mass eclipsing binaries with X-ray 
activity measurements from which we indirectly infer the Ha activity. Both samples yield consis- 
tent relations linking the amount by which an active object's temperature is suppressed, and its 
radius inflated, to the strength of its Ha emission. These relations are found to approximately 
preserve bolometric luminosity. We apply these relations to the peculiar brown-dwarf eclipsing 
binary 2M0535— 05, in which the active, higher-mass brown dwarf has a cooler temperature than 
its inactive, lower-mass companion. The relations correctly reproduce the observed temperatures 
and radii of 2M0535— 05 after accounting for the Ha emission; 2M0535— 05 would be in precise 
agreement with theoretical isochrones were it inactive. The relations that we present are appli- 
cable to brown dwarfs and low-mass stars with masses below 0.8 M© and for which the activity, 
as measured by the fractional Ha luminosity, is in the range —4.6 < log Lna/^boi ^ —3.3. We 
expect these relations to be most useful for correcting radius and mass estimates of low-mass 
stars and brown dwarfs over their active lifetimes (few Gyr) and when the ages or distances (and 
therefore luminosities) are unknown. We also discuss the implications of this work for improved 
determinations of young cluster initial mass functions. 

Subject headings: stars: low-mass, brown dwarfs — stars: fundamental parameters — stars: 
activity 



Introduction 



Observational evidence strongly indicates that the f undamental propert i es of low-mass stars can be 



altere d in the presence of strong magnetic activity (e.g. Morales et al. 2008 : Lopez-Morales! 2007 : Ribas 



2006). In particular, observations of numerous active, low- mass eclipsing binary (EB) stars have found the 
empirically measured stellar radii (R) to be inflated by ~10%, and the empirically measured stellar effective 
temperatures (T c g) to be suppressed by «5%, relative to the predictio ns of standard theoretical stellar 



evolution models, which better match the properties of inactive objects (see lCoughlin et al ho ill iKraus et al 
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2011 ; Morales et al. 201ot Torres et al. 2010 , and references therein). Because the mass-radius and mass- 
Teg relationships are central to our understanding of stellar evolution, resolving these discrepancies will be 



critic al to the ongoing development of accurate theoretical stellar models (for a discussion, see lStassun et al 



20101 ) . Accurate estimates of stellar radii are especially important in the context of searches for transiting 
exoplanets, which rely upon the assumed stellar radius/density to infer the planet radius/density. 

Activity effects also lead to errors in object masses (M) when these are derived from T e g. A particularly 
salient example is 2M0535 — 05, an EB in the Orion Nebula Cluster (age ~1 Myr) comprising two brown dwarfs 
(jstassun et al.ll2006i 120071 ). The primary and secondary brown dwarf (BD) components of 2M0535— 05 have 
dynamically measured masses of 6 0±3 and 39±2 Mj up , respectively, and T e ff ratio of T\jT-x — 0.952 ± 0.004 
(jGomez Maqueo Chew et al.ll2009n . That is, the system exhibits a reversal of the usual M-T e ff relation, such 
that the primary component is cooler than its companion. This behavior is not predicted by theoretical 
models for coeval BDs. Figure [T] s hows the 2M0 5 35— 5 system on the Hertzsprung-Russell (H-R) diagram 
compared to the 1 Myr isochrone of lBaraffe et al.l (|1998l ). The secondary BD's T e g and bolometric luminosity 
(J/bob calculated directly from the empirically measured T g and R) place it at a position that is consistent 
with that predicted by the model isochrone. In contrast, the primary is far displaced from its expected 
position, and so appears to have a mass of only ^ 2 5 Mj lln — more than a factor of 2 lower than its true 



mass — on the basis of its low T e g. iReiners et al.l ([20071 ) used spectrally resolved Ha measurements to 



show that, whereas the secondary BD in 2M0535— 05 is chromospherica lly quiet, the primary BD is highl y 
chromospherically active, perhaps a consequence of its rapid rotation ( Gomez Maqueo Chew et al. 2009). 
Thus, magnetic activity in the primary BD could be responsible for its highly suppressed T e g, similar to 
what has been seen for low-mass stellar EBs in the field. 

Given the exhibition of activity-related Toff suppression in 2M0535— 05, the first and only EB containing 
an active BD, it is likely that this phenomenon extends to other active BDs and low-mass stars in star forming 
regions and in the field. This has important implications for estimating masses from the H-R diagram, as is 
often necessary at the youngest ages. When the distance of a low-mass main-sequence star is known, mass 
can be estimated from its luminosity either via empirical mass-magnitude relations or model mass-Lboi 
relations, circumventing the need to use T e g. However, the use of luminosities can be problematic for very 
young objects, because they are sensitive to age, accretion, disk excess, and extinct ion — although careful 
evaluation of these parameters can mitigate these problems (e.g., iDa Rio et al.l I2012T ) . For low-mass stars, 
Tsff-derived parameters are useful on the pre-main-sequence, during which they evolve along nearly vertical 
tracks on the H-R diagram, and also while they remain active on the main-sequence lifetime. The activity 
timescale is relatively short at higher masses (<1 G yr for spectral typ es <M2) but increases substantially 
at lower masses (>7 Gyr for spectral types M5-M7) (|West et al]|2008l ). For substellar objects, determining 
physical properties is more complicated because such objects never reach a stable main sequence. Even with 
an accurately known distance and thus Tboi> the mass and age will be degenerate. Thus, mass estimates 
for field BDs are generally inaccessible in the absence of age information. At the youngest ages (<100 Myr) 
the mass-Lboi relationship for BDs is substantially flatter as they are still undergoing significant contraction 
of their radii, so the situation is more akin to pre-main-sequence stars where it is important to know T e g 
accurately. These evolutionary phases are also when BDs are most active, un like at field ages when they 



have reached spectral types >L4 for which activity is very rarely observed (e.g. Burgasser et al. 20021 ) . 



Since magnetic activity seems to alter the fundamental properties of both stars and BDs, it would be 
valuable to have an easily observable empirical metric with which to quantitatively assess the degree to which 
a given object's T e ff has been suppressed and its radius inflated. In the absence of a detailed understanding 
of the underlying physical causes of this effect, such a metric would be a useful "stop-gap" — it mitigates one 
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source of uncertainty in estimates of fundamental parameters for very low mass objects. The desired metric 
will allow us to compare magnetically active objects with their inactive counterparts and with non-magnetic 
evolutionary tracks in order to improve the mass estimates for young low-mass stars and BDs. 

The aim of this paper is to derive such an empirical metric by relating the degree of T e ff suppression 
and radius inflation to the strength of the Ha emission lin e, a commonly used and readily observable tracer 
of chromospheric activity ( Scholz et al. 2007 ; Berger 20061) . A challenge with any activity-based measure is 
that most active objects exhibit variability in their activity levels, and the amplitude of this variability can 
in some cases be quite large. For example, iBell et al.l (|2012l) found variability of up to 30% in Ha emission 
among a large sample of M0-M9 dwarfs in the field on timescales of minutes to weeks. Thus, the calibration 
of T c ff suppression and R inflation to Ha emission requires large statistical samples and/or calibration objects 
with highly accurate measurements in order to identify robust mean relations. Fortunately for this paper, 
which is concerned with the most highly active young objects which may experience the most significant T e a 
suppression and R inflation, stronger Ha emitters also tend to be less time- v ariable, with the typical Ha 
variability being less than 10% for the strongest Ha emitters ([Bell et al.ll2012l ). 



In Sec. [5] we describe our approach and the data we use to establish the empirical relationships of AT e g 
and AR vs. THa/^boi- In Sec.[3]we present the resulting relations and apply them to 2M0535— 05 as a test 
case, finding that the position of the chromospherically active primary BD in the H-R diagram can be fully 
explained as an offset from its theoretically expected position, due to the effects of activity on its R and T c s 
(Fig. [IJ. In Sec. 2] we discuss the broader application of these relations to other low-mass objects, and their 
possible ramifications for the inferred IMFs of young clusters. We conclude with a summary in Sec. [5] 



2. Methods and Data Used 

Our aim is to empirically determine a relationship between the temperature suppression (as compared 
to evolutionary models) and the level of activity as measured from Ha. We adopt a simple linear form for 
this relationship: 

AToff /T off = tot x (log Thq/ Tboi + 4) + b T (1) 

where AT e ff/T e ff is the fractional temperature offset (observed minus model), THa/Tboi is the ratio of 
luminosity in the Ha line to the bolometric luminosity, tot and 6t are linear coefficients, and the offset 
of +4 is for convenience given the typical value of log THa/Tboi ~ —4 for our study sample (see Sj3]). We 
determine a similar relationship for the fractional radius inflation: 

AR/R = m R x (log L Ha /L hol + 4) + b R . (2) 

In the following analysis, we use two different samples and approaches to determine the linear coefficients 
of Eqs. [Hand [2] In the first approach we use a large sample of stars without direct mass, radius, or T c ff 
measurements but with direct distances, spectral types, and reliable Ha measurements. The second approach 
uses the much smaller sample of stars in EBs that have directly measured masses, radii, and reliable T e ff, 
but for which we must use X-ray flux as a proxy for Ha. 



2.1. Nearby M Dwarfs with Ha Emission 



We first consider the large set of nearby field M dwarfs with well measured spectr al types and Ha 
equivalent widths (EWs) from the Palomar/Michigan State University catalog (PMSU; Reid et al. 19951 : 
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Hawlev et al. 19961) . In particular, following Morales et al. ( 2008 ). we restrict ourselves to the sample of 
746 stars with distances determined directly from trigonometric parallaxes. These distances allow Lboi to 
be calculated from the observed 2MASS Kg magnitudes together with if-band bolometric co rrections from 



Bessell et al 



(1998). The catalog reported spectral typ es are transformed to T e g according to lLeggett et al 



(|1996l ) and iBesselll (|l99lD . From these T eS and L hol ([Morales et all 120081 their Tables 1 and 2) we also 
calculate the stellar radii, R, according to the usual Stefan-Boltzmann relation. In order to isolate a consistent 
sample of main-sequ ence stars as represen tative as possible of Galactic disk population stars, we removed 
all stars identified bv iMorales et al.l (|2008l ) as being unresolved binaries, members of young moving groups, 
part of the old halo population, or evolved subdwarfs. The sample of 746 nearby stars with trigonometric 
distances is already dominated by field-age disk stars, so these cuts cause only 54 stars to be removed from the 
sample. We further restricted the sample to those with M < 1 M Q , leaving 669 stars in our sample. Finally, 



we use the solar- metallicity isochrone of Baraffe et al. ( 19981) with mixing-length a=1.0 to interpolate the 
predicted stellar masses (M), T c s , and R as functions of Lboi- Our results below do not change significantly 
for adopted isochrone ages in the range 1-5 Gyr, typical of field star ages for the Galactic disk population, 
so we adopt the isochrone at 3 Gyr throughout for simplicity. 

Fig- El shows the estimated T c g and R of the PMSU sample stars as a function of M. Strongly chromo- 
spherically active stars — defined here as those with Ha in emission — show a clear displacement to lower T e g 
relative to both the theoretical isochrone and to the non-active stars, whereas non-active stars more closely 
track the isochrone. The Ha active stars also show a displacement to larger R relative to both the theoretical 
isochrone and to the non-active stars. The mean offset of the Ha active stars relative to the isochrone (solid 
curve) is 10. Oct for R and —11. Id for T c ff, where a is the standard deviation of the mean (i.e., i.m.s./^/N). 
The mean offset relative to the non-active stars (dashed curve) is 6.8<r for R and — 7.5<t for T c g. The 
non-active stars themselves show a mild displacement to lower T e g and larger R relative to the theoretical 
isochrone. We note that at least some of the "non-active" stars in the sample ma y in fact possess mild 
chromospheric activity . Mildly active M-type stars can exhibit Ha in absorption (e.g. IWalkowicz fc Hawlev 



2009: 



West et al.l 1201 II ) and thus would not be identified as strongly active according to our Ha emission 



criterion but could still manifest mild activity-related effects. Alternatively, these offsets for the non-active 
stars may simply indicate mild systematics in the transformation from the observed if -band magnitudes 
into masses and radii via bolometric corrections; indeed, the polynomial fits to the non-active stars closely 
parallel the theoretical isochrone (except at the highest masses; see also Sec. I3.1j) . In what follows we con- 
servatively measure the T e g and R offsets of the Ha active stars relative to the non-active stars, effectively 
using the latter to calibrate the derived stellar masses and radii to the theoretical isochrone. For reference, 
the polynomial fits with which we describe the non-active stars (dashed curves in Fig. [5]) are: i?/R© = 
X^o^M/M©) 1 and T cff /K = Y,t= U(M /M Q y , with n = {1.3836,-10.9203,37.7554,-52.3712,26.8055} 
and U = {-1384.1, 41857.5, -136147, 192515, -96288.4}, applicable for M < 0.8 M Q . 



Fig. 3(a) shows the offsets in estimated T c g and R for the Ha active sample stars as a function of the 
observed EW(Ha). We fit a simple least-squares linear relation to each (dashed lines), yielding the following: 

AR/R [%] = (1.4 ± 0.8) x EW(Ha)[A] + (6.1 ± 3.2) 

AT cff /T cff [%] = (-0.6 ±0.3) x EW(Ha)[A] + (-3.0 ± 1.3). 

For example, for EW(Ha)=4A, the relation predicts a T c g suppression of «5% and a radius inflation of sall%. 
There is substantial scatter in the data, nonetheless a Kendall's r correlation test gives a null-hypothesis 
probability of only 0.9% and 0.7% that the quantities are not correlated in AR and in AT c g, respectively, 
indicating that the correlations are statistically significant at >99% confidence. For comparison, a linear fit 
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to the T e fj and R offsets cafculated refative to the theoretical isochrone is shown as a solid line. The linear 
fit cofhcients are nearly identical to those above, except of course for a larger mean offset in both quantities. 

To calculate the Ha luminosity, Lna, from the EW(Ha) measurements for Equations [THU we determined 
scaling factors from a set of stellar atmosphere models for T c g ranging from 3000 to 5000 K in steps of 100 K 
and log g = 4.0. We use the STARd ustyl999 model spe ctra, which are based on the NextGen models updated 
with new H 2 and TiO opacities ([Allard et al.l 12000). We obtain the continuum flux around the Ha line 
(6500-6600 A) for each T c g, and we fit the flux-T^ relation with a third-order polynomial. The measured 
EW(Ha) is then multiplied by the Ha continuum flux from this relation, and multiplied by the surface area 
of the star using the measured R calculated as described above, thus giving the total Lhq- Fig. |3(b)| shows 
the same data as in Fig. 3(a) but now with EW(Ha) converted to I/Ha/iboi- The resulting relations between 



the T c e and R offsets vs. Ln a /Lhoi are discussed in Sec. 



2.2. Low-mass Eclipsing Binary Stars with X-ray Emission 



In the second approach, we use the small se t of lo w-mass EBs with accurately measured M, R, T e ff, 
and X-ray luminosities (Lx) from Lopez-Moralesl ( 20071) Fl The sample includes 11 individual stars in 7 EB 
systems spanning the range T o ff=3125-5300 K and M =0.21-0. 96 M Q . 

We begin with the correlation of AR vs. Lx / £t>oi alr eady demonstrated in that work, which we rederived 
(F ig.|4(ajj) using the fu ndamental stellar data compiled i n lLopez-MoralesI (|2007T) and the same 3 Gyr isochrone 
of Baraffe et al. ( 19981 ) as above. Lopez-Morales ( 2007 ) did not discuss the complementary correlation with 
AT c fj, but this information is also contained in the EB d ata, so in Fi g.|4(a)|we also derive the the relationship 
A T R ff vs. Lx /£hni, again using the data compiled in lLopez-MoralesI ([20071 ) and the 3 Gyr isochrone of 



Baraffe et al.l (|1998l ). Linear fits are shown for both relations: 



AR/R [%} = (15.5 ±2.9) x log L x /L bo i + (57 ± 9) 
AT cff /T off [%] = (-6.2 ± 3.2) x log L x /L hol + (-23 ± 10). 

A Kendall's r correlation test gives a null-hypothesis prob ability of 0.2% for th e correlation with AR, 
indicating a significant correlation at 99.8% confidence (see also Lopez-Moralesl 2007 ). The correlation with 
AT g is not as strong, but is nonetheless modestly significant at 94% confidence according to the Kendall's 
r test. Interestingly, the AR and AT e ff relations, which for the EB sample are measured independently, very 
nearly offset one another in terms of their effect on Lboi- For example, at log Lx/L bo \ = —3, the relations 
give AR ss 11% and AT e ff ss —5%, implying ALboi ~ 1%- Thus £t>oi is an approximately conserved quantity. 

To convert the observed LxjL\, \ into Ln a /L bo i for use in Equations [TH2j we use published activity 
data for two different sam ples having both X -ray and Ha emission measurem ents (see Fig. [5]) : low-mass 



stars in young associations ( Scholz et al. 2007 ) and active M dwarfs in the field ( Delfosse et al. 19981 ). Fig. [5] 



demonstrates that there is a robust correlation between these two quantities over a wide range of activity 
levels, albeit with considerable scatter. A linear least-squares fit to this correlation is overplotted as solid and 
da shed line for the first an d second sample, respectively. The activity data for M-type field dwarfs published 
bv iReiners fc Basril (|2007l ) shows a similar trend (circles in Fig. [5]). We note that these correlations can 



1 We use the "case 1" Lx values from lLopez-Mor alcs (2007). The results do not change significantly if we adopt the "case 2" 
or "case 3" Lx values instead. 



- 6 - 



be extended to lower activity le vels, as they are consistent with the simultaneous activity data for BDs 
published bv lBerger et al.l (|2010l . their Table 5). For the young low-mass stars sample, the fit is: 

\og{L x /L hol ) = (1.5 ± 0.7) x logLna/iboi + (3.0 ± 2.4). 

For the field M dwarfs and BDs, the fit is: 

log(L x /L hol ) = (1.1 ± 0.3) x logWLbd + (1.0 ± 1.1). 

For our analysis we have cho sen to use the second relation because it better reflects the range of X-ray 
activity levels observed in the Lopez-Moralesl ( 2007 ) sample of low- mass EBs (see Fig. 4(a)). This is also 
a more conservative choice, in that it associates a lower Lx/L\, \ for a given Ln a / Lt, i, and therefore will 
predict smaller absolute AT g and AR offsets for a given observed Ln a /Lb \. 

Fig. |4(b)| shows the same data as in Fig. 4(a) but now with Lx/L\, Q \ converted to Ln a /Lboi- The 
resulting relations between the T c g and R offsets vs. ina/^bol are discussed in Sec. [3] 



3. Results 

3.1. Empirical Relations Linking T e g Suppression and Radius Inflation to Ha Emission 

We have determined empirical relationships linking the degree of temperature suppression, AT e s/T e s, 
and radius inflation, AR/R, to the fractional Ha luminosity, LHa/^boi, using two independent samples and 
methods (see Figs. |3(b"J| and [4(b)] ). The fit coefficients from Eqs. [TH2]for the two samples are summarized in 
Table [TJ The fit coefficients are defined such that the left-hand side (LHS) of Eqs. [TJ[2]are in percent units. 

The AT e g and AR fit coefficients in Table Q] are consistent to within l-2cr between the two samples, 
although the formal uncertainties on the fit parameters are large for the field M-dwarf sample. This likely 
reflects the large scatter in the Ha measurements for that sample, perhaps stemming from intrinsic stellar 
variability. In addition, for the field M-dwarf sample the determination of Ln a /L\, i involved several cal- 
culated quantities (i.e., i?, Lboi) whereas for the EB sample these stellar parameters are measured directly 
and accurately. Indeed, for the field M-dwarfs the statistical significance of the AT e g and AR correlations 
was stronger before we converted the directly measured EW(Ha) to I/Ha/^boi- Even so, both the AT c g and 
the AR relations for both samples are confirmed to be significantly correlated with Ln a /Lb a i according to 
a Kendall's r test (Table [TJ). Being a non-parametric rank-correlation test, Kendall's r does not depend on 
the assumed functional form of the relationship, and thus robustly indicates the presence of a correlation 
even if the significance of the assumed functional parameters is modest. 

For our final best-fit relation, we calculated the weighted average of the fit parameters for the two 
samples. The resulting best- fit coefficients for Eqs. [TJ and [5] are also listed in Table [TJ and the corresponding 
final fits shown in Figure [6] As is evident in Fig.[6j the r.m.s. scatter of the data about the mean relations is 
large (4.2% and 9.2% for AT e g and AR, respectively), driven primarily by the scatter in the field M-dwarf 
data (see Sec. l2.1|) . Nonetheless, the statistical significance of the final combined relations is strong (Table[lJ. 

The input samples used to determine these relations include low-mass stars with 0.2 < M < 1 M Q . 
There is no evidence that the T e g suppression and R inflation effects change qualitatively at lower masses 
(see Fig. [2]), and we are able to validate them successfully at BD masses (see below). However, the relations 
appear to diverge at masses >0.8 M (see Fig. [2]) and so we caution against their use at such high masses. 
Given the range of Ln a /L\, i spanned by the input data samples, these relations are applicable for —4.6 < 
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Table 1. Fit parameters for Equations Q] and [2] 





m b 


t conf. a 




Field M-dwarfs b 


ATeff 
AR 


-3.12 ±3.15 -5.1 ±0.7 
8.00 ±7.63 11.2 ±1.6 


92.9% 
91.2% 


Eclipsing Binaries 


ATeff 

AR 


-6.64 ±3.47 -3.0 ±1.0 
16.64 ±3.15 6.6 ±0.6 


94.3% 
99.8% 


Averaged Final Relation d 


ATeff 

AR 


-4.71 ±2.33 -4.4 ±0.6 
15.37 ±2.91 7.1 ±0.6 


96.0% 
98.6% 



Statistical confidence of correlation from 
Kendall's r test. 

b See Fig. [3(b)] 

c See Fig. [4(b)] 

d See Fig. El 

Note. — The fit parameters are defined such 
that the LHS of Equations [T]-[5] are in percent 
units. 
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log Lua/L^oi < —3.4 (see Fig. [6]) , the upper limit corresponding appr oximately to the empirical "saturation " 
level of log ina/iboi ~ —3.3 observed in low-mass stars and BDs (e.g. Barrado v Navascues fc Martinll2003h . 
Therefore the relations likely cannot be extrapolated to log iHa/iboi > —3.3. Note also that for log I/Ha/^boi 
< —4.6, Eqs.[T][2]give positive AT e ff and negative AR, which is likely not physical. More likely the AT ff and 
AR offsets simply approach zero at very low activity levels, and thus we caution that our relations should 
not be extrapolated to arbitrarily low Ln a . 



3.2. Application to 2M0535 05 



2M0535— 05 is to date the only known system of BDs with directly and accurately measured masses, radii 
T e g, and L^ a , all at a well constrained system age ( Stassun et al. 20061 2007 ; Gomez Maqueo Chew et al 



2009). 2M0535— 05 is therefore an important empirical test case for assessing the efficacy of the activity-based 
AT c ff and AR relations that we have determined above. The known masses of the primary and secondary 
BDs in 2M0535-05 are 60±3 and 39±2 M Jup , respectively. 

To demonstrate how the apparent masses of the 2M0535— 05 primary and secondary BDs are altered 
by our empirical relations, we calculate how the observed positions of the two BDs in the H-R diagram 
are altered by adjusting the observed T e g and R using Eqs. [TH2]and Table [T] In effect, we are using our 
empirical relations to show how the 2M0535— 05 system would appear in the H-R dia gram were the s ystem 
completely inactive. The spectral type of M6.5±0.5 for 2M0535— 05 determined by Stassun et al. ( 20061) 
from high-resolution ff-band spectroscopy i mplies an average T e g=2688±55 K, based on the near-infrared 
spectral- type-Tgff relation of ISlesnick et al. ( 2004). Weight ing this average T e g by the i?-band primary-to- 



of Ti/T 2 = 0.952 ± 0.004 (|G6mez Maqueo Chew et al 



secondary flux ratio of 1. 6:1 found bvptassun et al. ( 2007 ) and using the accurately determined T e s ratio 

; ff =2640±60 and T cff =2770±60 K 
oge 

2009) then give the individual 



2009) then gives T ( 



for the primary and second ary components, respectively. These individual T ? e tog ether with the accurately 
measured individual radii (jStassun et al.ll2007t iGomez Maqueo Chew et al. 
iboi- These observed T e g and -Lboi for the primary and secondary components are represented in Fig. Q] 
as blue symbols, from which one would infer masses (using the Baraffc ct al. evolutionary models) of «27 
and «40 Mj up , respectively, based on the observed T ff. That is, the inactive secondary is inferred to have 
approximately its true mass, but the active primary appears to have a mass that is a factor of 2 lower than its 
true mass, and moreover the primary appears to be much younger than its (presumably coeval) companion. 



We use the Ln a /L^ i measurements of iReiners et al.l (|2007l) . who found log Ln a /L^ i — —3.47 for the 
active primary and log iHa/^boi < —4.30 (upper limit) for the inactive secondary. From Eqs. [1] and[2J we 
find that the primary has been displaced from its inactive T e g by an amount AT c ff = —6.9 ± 1.4%, and 
displaced from its inactive R by an amount AR = 15.2 ± 1.7%. For the secondary, the resulting AT e ff is at 
most —3.0 ± 0.9%, and AR is at most 2.5 ± 1.1%. For both BDs the AT D s and AR displacements nearly 
cancel such that Lboi is preserved to within 1-3%. 

As shown in Fig. [TJ shifting the H-R diagram position of the primary BD according to its calculated 
AT c ff and AR shows that, if it were inactive, it would be in remarkably good agreement with the theoretically 
expected position for its known mass. The observed Ln a /Lb \ of the active primary is near the "saturation" 
limit and therefore near the upper limit of observed act ivity levels in yo ung low-mass objects, for which 
the Ha emission is found to be mostly non- variable (e.g. iBell et al.l 120121 ). Thus we do not expect a large 
contribution to the uncertainty of the active primary in Fig. [T] due to Ha variability. We do not shift the 
position of the secondary as its Lhq is an upper limit only. 
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Discussion 



4.1. Ha, magnetic activity, and cooled/inflated stars 

In developing these empirical relationships, we posit a correlation between chr omospheric activity, ra dius 
inflation, and T B ff suppressio n, where activity is implicitly tied to magnetic fields. Chabrier et al.l |2007h and 



MacDonald &: Mullanl ([20091 ) have suggested that a sufficiently strong field could suppress convection, inhibit 
heat transfer, and thus inflate (and cool) the stellar surface. Since such a field would also likely result in 
chromospheric activity, one might therefore expect the correlations that we have derived. However, there are 
two significant outstanding questions in creating the above link: (1) Is the relation between Ha emission, field 
strength, and radius inflation / T e g suppression monotonic? (2) Are the magnetic fields in fully convective 
objects strong enough to account for the radius inflation and T c ff suppression? 

First, it remains unclear how well Ha emission correlates with the strength of the magnetic field. 
Reiners fc Basril (|2007T) have attempted to relate field strength and Ln a /L hol using the FeH line. They find 
a correlation between the two, but with different power law scalings for different spectral types. Whether this 
is an age or mas s effect remains uncert ain. While this correlation is encouraging for the empirical relations 
we provide here. iReiners fe Basril (|2007j) use slowly rotating field stars, whose magnetic properties may differ 
from the most active (younger) objects in which we are interested. As noted previously, Ha em ission is also 
variable, though the most active objects typically have the lowest variability (jBell et al.ll2012l ). Concurrent 
multi-wavelength monitoring of multiple chromospheric trace rs including radio and X-ray, sh ow that the 

Thus it is 



individual tracers do not always vary on the same timescales (jBerger 



200t 



(J Berger 



et al 



200 



unclear that any individual wavelength range is representative of the total energy contained in the field. 
Despite these different trends in variability, that Lx also correlates with T e g suppression and radius inflation 



(see Fig. 4(a) ) corroborates our use of chromospheric activity to account for the displacement of objects from 
model isochrones (see also, e.g.. Morales et al.ll2008l and references therein). 



The second concern, whether or not magnetic fields are sufficiently powerful to inflate the star s and sup- 
press their surface temperatures, is currently impossible to address without appealing to models. iBrownine 
([20081 ) has shown in global numerical models that fully convective stars can host large Kilogauss strength 
fields, but t hese alone are too weak to p roduce the observed r adius inflation and T c ff suppression in 
2M0535-05 (jMacDonald & Mullanl l2009h . Ichabrier et~all (|2007h also suggest that fully convective ob- 



jects should be less affected by the same convective inefficiencies invoked to explain radius discrepancies 
at higher masses. It is possible that a combination of rotation and magnetic activity contribute to both 
inflation/suppression and Ha emission in such a way as to produce our empirical relation without a causal 
correspondence between Ha and the magnetic field. 



An alternative expl anation for star s with T e g deficits is a spot covered surface (e.g. lLacvlll977j) . For 



the case of 2M0535— 05. iMohantv et al.l (|2010h have argued that a model for the active primary with 70% 



(axisymmetric) spot coverage could explain the peculiar mass-temperature relationship in that system. Since 
spot coverage is also controlled by magnetic fields, a correlation with Ha might still exist. It remains unclear 
how to interpret spot coverage o f greater than 50%; perhaps the analogy with solar type spots breaks down 
in such extreme systems. Indeed. IMohantv fc Stassunl (|2012l ) have now shown from a spectral fine-analysis of 
the 2M0535— 05 system, observed at high resolving power during both primary and secondary eclipses, that 
such a large-spot scenario is strongly disfavored as an explanation for the T c g suppression of the primary BD 
in 2M0535— 05. Thus, while the T e g suppression mechanism produces a clear correlation with chromospheric 
Ha activity as we have shown here, it evidently does not in all cases effect this correlation directly through 
surface spots. 
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4.2. Impact on estimates of object masses 



The vast majority of masses for stars and BDs in young clusters can only be determined by com- 
parison with theoretical evolutionary tracks in H-R diagrams, using either luminosity or T e g or bot h (e. 



Kenvon fc Hartmannlll995t iHillenbrandl 1 1 997t iLuhman et al.l ll998 



Moraux et al. 2007; Scholz et al. 201 



Whi le direct mass measurements are available for calibration of these models at higher masses (M > 0.3 M Q ; 



e.g., iHillenbrand fc White! 120041 ) . at the lowest stellar masses and into the substellar regime these models 
are essentially untested. Therefore, if there is indeed a relation between magnetic activity and R inflation / 
T e g suppression, this will also affect estimates of stellar and substellar masses derived from T e ff, especially at 
young ages when activity levels are high. In the following, we examine two ways of estimating masses that 
are commonly used in the literature and investigate the impact of magnetic activity on the derived masses. 

First we use T e g to estimate masses from model isochrones at 1 Myr. At each model mass, we consider 
a range of activity levels and apply offsets to the model T e g based on our Eq. [T] (using "averaged final" 
coefficients from Table [lj . We then use these suppressed T c g values to estimate the masses that would be 
inferred from the isochrone. The resulting curves are shown in the upper panel of Fig. [7J (Note that for 
log inc/iboi = —4.7 the offset with Eq. Q] becomes zero .) We consider mas ses of 0.02-0 .60 Mq (i.e., T^ff 



2500-3500 K), for which we use a combination of BCAH (jBaraffe et al.lll998l) and Dusty (jAllard et al.lboOll ) 
isochrones. (Using BCAH models at masses >0.1 M and Dusty models at lower masses, we fit a third order 
polynomial to mass as a function of T e g.) 

As expected, this procedure leads to a systematic underestimation of the masses. At high levels of 
activity, the effect can be substantial. For log Ln a /Li, i = —3.3 , which corresponds to the saturation limit 
in young low-mass stars (IBarrado v Navascues fc Martini 120031) and the upper limit of the observed Ha 
activity in young associations (jScholz et al.l 120071 ). the mass estimates are a factor of ~ 2 lower than for 
objects with low levels of magnetic activity (log L-a a /L^ \ < —4.5). In Fig. [7]we overplot the datapoint for 
the primary component of the eclipsing BD 2M0535— 05. As discussed above, its mass would be estimated to 
be 0.03 M Q , whereas its true mass is about twice as high. We note that because the mass-T e ff relationship 
is less steep at older ages, the underestimation of mass is less severe at older ages (>100 Myr). 

In the second test case, we use Lboi derived from i^-band absolute magnitudes to derive masses. In this 
case the influence of magnetic activity is introduced by the T e ff dependence of the bolometric correction, since 
the bolometric luminosities are (as outlined above) practically not affected. We consider the same range of 
m asses as in the fir st case, and we compute spectral types from the suppressed values of T c ff usi ng the relation 
of"" 



Luhmanl ([2003, fitted linearly) . We then computed -fT-band bolometric corrections from ISlesnick et al 



(200 



3) from the spectral types, and we combine these with the model K-band absolute magnitudes to find 



Lboi (we used a linear fit of model Mr as a function of T c g). Finally, we estimated masses from the model 
isochrone and £t>oi, and the result is shown in the lower panel of Fig. [7] 

This second method still leads to an underestimate of the masses, but the effect is much smaller than 
when directly estimating masses from T c g. The curves are almost flat, i.e. the estimated masses depend 
insignificantly on the level of magnetic activity. The change of the bolometric corrections with increasing 
magnetic activity is quite small (at most 10%), resulting in relatively minor changes of a few percent in the 
mass estimate. With this method, the mass for the primary component of the eclipsing BD 2M0535— 05 
would be 0.054 M and thus only marginally smaller than the true value. 



These effects are systematic biases in estimated masses due to T e g suppression, so this will always be 
significant for statistical studies of populations (e.g., in young clusters). For individual objects, this bias 
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may be somewhat lower than the uncertainty in T e g due simply to our limited ability to accurately model 
spectra of 2500-3500 K objects. For example, the best benchmark objects currently available over the 
2500-3000 K range — very low-mass stars and BDs with dynamical mass me asurements — show that model 
atmospheres likely harbor large systematic errors («250 K; IDuduv et "alboioh . This uncertainty corresponds 



to a fractional error in T e g of ~8%-10%, which is somewhat larger than the typical offset we find due to T e s 
suppression. 

An important application of T c g-based mass estimates is the determination of mass functions (MFs) for 
young clusters. Typical early to mi d M-type dwarfs with ages from 10 to 100 Myr exhibit Ha emission in the 
range log L Ha /L hol = -4.2 to -3.3 (jBarrado v Navascues fc Marti'nll2003l:IScholz et al.ll2007h . For BDs at the 



same ages the available data are sparse but indicate an upper limit around —3.7 (jBarrado v Navascues fc Martin 



20031 ). Thus, the upper panel of Fig. [7] implies that objects inferred to be BDs from T e ff-based (or spectral 
type based) mass estimates could actually be low-mass stars since masses will be underestimated by up to 
a factor of 2. 

We further evaluated the impact of this effect on measurements of the slope of the mass function, 
a (in dN/dM cx M~ a ). For this purpose, we assume a measured slope of a = 0.6, which is consistent 
with a number of studies in very young clusters (see Scholz et al. 2012, submitted, and references therein), 
calculate the mass function based on that slope for M < 0.6 M Q , correct the masses for a given level of 
Ha emission using the results given above, and re-determine the slope a. If magnetic activity and thus 
the level of Ha emission is constant across the low-mass regime, a is practically unchanged, because all 
masses will be underestimated by about the same factor. Based on the available measurements, however, 
i t seems more realistic to assume tha t BDs have on average a lower level of activity than low-mass stars 



(jBarrado v Navascues fc Martrnll2003l ). If we start with log Lna/^boi = —4.0 for BDs and —3.5 for low-mass 
stars, the masses have to be corrected by factors of 1.3-1.7 for BDs, and by 2.2-2.5 for low-mass stars. As 
a result, the slope changes from a = 0.6 to a = 0.5. A more extreme mass dependence of the activity level 
will enhance this effect. In general, we expect that a will be underestimated by > 0.1, if the masses are 
estimated from T e g and activity is not taken into account. In addition, the peak of the mass function could 
be underestimated by up to a factor of 2. For a more quantitative assessment of these effects a larger sample 
of Ha measurements for young stars and BDs is needed. 



5. Conclusions 

In this paper, we have shown that there exists a correlation between the strength of Ha emission in 
active M-dwarfs, and the degree to which their temperatures are suppressed and radii inflated compared with 
inactive objects and theoretical evolutionary models. By applying these relations, we are able to infer the 
amount by which an active objects' temperatures have been suppressed, and thereby improve the accuracy 
of estimates for their masses and radii. We use the brown dwarf eclipsing binary system 2M0535— 05 as a 
benchmark for our model. We expect these relations to be most useful for corr ecting estimated masses and 
radii of low-mass stars and brown dwarfs over their active lifetimes (few Gyr; West et al. 20081) and when 



the ages or distances (and therefore the luminosities) are unknown. We have shown that failing to account 
for the effects of activity can cause significant errors in estimates of stellar and substellar masses derived 
from T c ff, and smaller, but systematically biased errors in temperature and radius. 

If these empirical corrections are corroborated by future observations, they will prove valuable not 
only for individual objects, but also for studies of stellar populations. In the case of individual objects, 
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reliable stellar properties are invaluable for exoplanetary studies, where exquisite knowledge of the host star 
is required to infer planet properties. For very young stellar populations, where activity levels are highest, 
we have shown how underestimated masses can substantially shift the inferred initial mass function. Such a 
change would necessitate revisions to star formation models and population synthesis models because, e.g., 
the observed fraction of brown dwarfs and low-mass stars might be substantially altered. 

While promising, the correlations we have derived contain significant scatter, and they are currently 
limited by the lack of a single sample of stars with both Ha and direct radius measurements. A sample of 
Ha measurements for objects with directly measured radii and temperatures will allow us to better assess 
our relations and determine if the scatter is intrinsic or if it is caused by the intermediate steps necessary 
in constructing our relation (e.g., converting Lx to Ln a ). Despite this limitation, we have chosen to pursue 
a correlation with Ha emission (rather than X-ray, or radio, for example) because of the relative ease of its 
measurement even in the substellar regime. In principle, a reliable activity tracer in the near-IR would prove 
even mor e useful by making m easurements easier for cooler objects. No such tracer has yet b een identified 



(see, e.g., ISchmidt et al.l 120121) . although the He I line at 10830A may be a possibility (e.g. iDupree et al 
1992h . 



In order to make progress on understanding how T e g suppression and radius inflation relate to chromo- 
spheric activity, a larger sample of EBs with precise masses, radii, and Ha measurements are needed. Toward 
this goal, we encourage other researchers to publish Ha measurements for their targets, as this is usually 
readily available from the spectra used to determine radial velocities. Ind eed, a number of low-mass EBs 
with accurate masses and radii have been published in the last few years (e.g. 



Vaccaro et al 



2007 



Irwin et al 



Irwin et al 



2009[ iMorales et aUboid iKraus et al1l201l[ iHelminiak et alJboilt iHelminiak fc Konackill2011 
2011), potentially increasing by a factor of 2-3 the small sample that we have used from Lopez-Morales 
([20071 ) . We are currently collecting new Ha measurements for these published EBs that lack resolved Ha 



measurements in order to improve the empirical relations that we have presented here. 

Finally, the relations we have determined already indicate quite clearly that the radius inflation and 
temperature suppression mechanism operates in such a way that the temperature suppression and radius 
inflation almost exactly cancel in terms of their effect on the bolometric luminosity. Moreover, the relations 
between activity, T e s suppression, and radius i nflation do not appe ar to manifest any obvious discontinuity 
across the fully convective transition (see also lStassun et al.ll2010i ): the followup observations we have un- 
derway should help to refine this. These are important, fundamental clues to the physical nature of these 
effec ts, and should help to constrain theoretical models that are being developed to explain these phenomena 
(e.g. IChabrier et~aI1l2007t iMacDonald fc Mullanlliooj ). 
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Fig. 1. — Hertszprung- Russell diagram for the primary and secondary c omponents of the b rown- 



2006, 2007 



dwarf eclipsing binary 2M0535— 05 in the ~1-Myr Orion Nebula Cluster (St assun et al 
Gomez Maaueo Chew et al.l 120091 ) . The measured T e g and Lboi (the latter calculated from the directly 
measured T e g's and radii; see Sec. l3.2p for both brown dwarfs are represented as blue symbols. Measurement 



uncertainties in T c g and Lboi are represented by the error bars at upper right. The dynamically measured 
masses of the primary and second ary are 60±3 and 3 9±2 Mj up , respectively, represented as red bars on 
the 1 -Myr theoreti c al iso chrone of Baraffe et al. ( 19981 ). The measured Lna/iboi for the two components 
from lReiners et alj (|2007l ) are indicated next to the blue symbols. The inactive secondary appears close to 
its expected position on the isochrone, whereas the active primary appears far cooler than expected. The 
primary therefore appears to be much younger than the secondary and to have a mass of only «27 Mj up 
based on its observed T G ffj a factor of 2 lower than its true mass. Shifting the position of the active primary 
(arrow) using our empirically calibrated Ha-based relations for T e g suppression and radius inflation brings 
the primary into much closer agreement with its theoretically expected position in the HR diagram (black 
symbol); this is where the active primary would be if it were not active. This figure is shown in color in the 
electronic version only. 
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Fig. 2. — Effective Temperature vs. Mass (left) and Radius vs. Mass (right) for M-dwarfs with trigonometric 
distances and Ha measurements from the PMSU catalog. Active objects (i.e., thos e with Ha in e missio n) 
are represented by filled symbols. For reference, the 3 Gyr theoretical isochrone of iBaraffe et al.l (|1998T ) is 
represented in both figures as a solid curve. A polynomial fit to the non-active objects is represented by a 
dashed curve. The Ha-active dwarfs are found to be significantly displaced to lower T e ff and larger radii as 
compared to both the theoretical isochrone and the non-active dwarfs. 
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Fig. 3. — Suppression of effective temperatures (left) and inflation of radii (right) as a function of Ha emission 
strength for the Ha-active sample from Fig. [2] (filled symbols), measured relative to the non-active stars. In 
(top), the abscissa is the directly measured Ha equivalent width, in (bottom) the Ha equivalent widths have 
been converted to fractional Ha luminosity. In all panels, the solid line is a linear fit to the T e g and radius 
differences relative to the theoretical isochrone (solid curves in Fig. [2J ; the dashed line is a linear fit to the 
T e s and radius differences relative to the non-active stars (dashed curves in Fig. [5]). For our analysis we use 
the more conservative T e g and radius differences measured relative to the non-active stars (shown here as 
filled symbols and dashed lines). 
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Fig. 4. — Su ppression of effective t emperatures (left) and inflation of radii (right) for low-mass eclipsing 
binaries from Lopez-Moralesl ( 2007 ) as a function of fractional X-ray luminosity (a) and as a function of 



fractio nal Ha luminosity (b). The AT e g and AR are measured relative to the 3 Gyr isochrone of lBaraffe et al 
(|1998l ). The fractional Ha luminosities are based on the fractional X-ray luminosities, using the empirical 



X-ray-to-Ha luminosity relation in Fig. [5] In all panels, the dashed line is a linear fit to the data. 
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Fig. 5. — Fractional X-ray luminosity vs. fraction al Ha lumin o sity f or low-mass star samples from field 



dwarfs ( Dclfosseetal 



19981 crosses), young stars Scholz et al. ( 2007 . plus symbols), and field M-dwarfs 



(jReiners fc Basrill2007L circles). Linear fits to the low- mass star samples are shown as solid (fit to plus 



symbols) and dashed (fit to crosses) lines. 




Fig. 6. — Temperature suppression (left) and radius inflation (right) as a function of fractional Ha luminosity 
for both field M-dwarfs (filled symbols) and low-mass eclipsing binaries (error bars). Relations from Figs. |3(b)j 
and |4(b)| are reproduced (dashed lines). The final averaged best-fit relation in each panel is shown as a solid 
line. See Table [T] for the linear fit coefficients to these relations (Eqs.[T]and[2|). 
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Fig. 7. — Inferred masses of low-mass stars and brown dwarfs as a function of Ln a /Lb \, derived using the 
BCAH98 and DUSTY evolutionary tracks and the empirical relation between T e g suppression and magnetic 
activity derived in this paper. In each panel, the true object masses are labelled on the curves, whereas 
the masses that would be inferred at a given activity level are shown on the ordinate. At low activity 
levels (log Lua/Lboi ^ — 4.6), the inferred mass is close to the true mass. Top panel: Masses are inferred 
directly from the observed T e g, which is susceptible to suppression at high activity levels, and therefore the 
inferred masses are systematically lower than the true masses. Bottom panel: Masses are inferred from the 
bolometric luminosity, estimated from the if-band absolute magnitude and bolometric correction (as would 
be appropriate for an object of known distance). Here the T c g suppression effect enters only weakly via 
the T e ff-dependency of the if -band bolometric correction. The dash-dotted line in each panel shows the 
substcllar limit. The triangle in each panel shows the position of the primary component in the eclipsing 
binary brown dwarf 2M0535— 05, which has a true mass of 0.06 M©. 



